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X. .4 new Method of deducing a first Approximation to the Orbit 
of a Comet from three Geocentric Observations. By James Ivory, 
A. M. Communicated by Henry Brougham, Esq. F.RS. 

Read February 17, 1814. 

Comets are distinguished from the planets not only by the 
peculiarities that immediately strike the eye, but likewise by 
the circumstances attending their motion in the heavens. All 
the planets move round the sun in orbits nearly circular; they 
never deviate far from the ecliptic on either side ; and they 
move in a manner not extremely irregular, and in one direc- 
tion, according to the order of the signs in the zodiac. Comets, 
on the contrary, when they first come into view, assume gra- 
dually greater degrees of brightness, which they again lose 
by like gradations, and then disappear ; thus seeming to visit 
the neighbourhood of the sun for a short time only, after which 
they retire into the immensity of space : they are seen in all 
quarters of the heavens : and their motion is exceedingly va- 
rious and irregular ; confined to no direction ; sometimes 
greatly curved, and often hardly distinguishable from a rec- 
tilineal course. If, to phenomena so dissimilar, we add the 
prejudice which almost universally prevailed, that comets 
have only a temporary existence, and are produced by occa- 
sional causes, we shall not perhaps have much reason to be 
surprised that the true account of those bodies, which repre- 
sents them as forming a part of the same system with the 
planets, eluded the sagacity of Kepler, to whom we are 
mdcccxiv. R 
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indebted for the first accurate knowledge of the laws of the 
planetary motions. This step in our knowledge of the uni- 
verse was reserved for Sir Isaac Newton. The principal and 
leading discovery of that great philosopher consisted in gene- 
ralizing the laws of Kepler ; in proving that they are neces- 
sary consequences of a more general fact, namely, that all the 
planets are continually deflected from a rectilineal motion to- 
wards the sun in the inverse proportion of the squares of their 
distances from that body. He demonstrated that the motions of 
such a system of bodies must be performed in the conic sections, 
having the sun in the focus, the species of the curve depend-, 
ing upon the proportion of the rectilineal velocity to the quan- 
tity of the deflection towards the common centre. This theory 
comprehends an infinite variety of motions, all flowing from 
one common principle ; and the ellipse alone, by the changes 
of form which it undergoes according to the degrees of its 
eccentricity, seems, at one extreme, when it is greatly elon- 
gated, as well adapted to account for the phenomena of the 
comets, as it is, at the other extreme, when it differs little from 
a circle, to represent the motions of the planets. 

To try if this theory will account for the actual appearances, 
it is necessary to determine, by means of terrestrial observa- 
tions, the magnitude and position of the curve in which a comet 
moves round the sun. Having selected as many observed 
places of a comet as are necessary for this purpose, the re- 
maining observations will serve as so many tests of the accuracy 
of the theory. When a comet appears for the first time, it is, 
indeed, hardly possible to determine its orbit with exactness. 
The very eccentric ellipses in which those bodies move, allow 
them to be seen only when they are near the sun ; and in this 
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situation, the nicest observations do not enable us to distinguish 
the real orbit from the parabola with which it intimately co- 
incides. It is only after one or more reappearances of the 
same comet, that we can hope to discover the period of its 
revolution round the sun, and the mean distance of the ellipse 
in which its motions are really performed : unless indeed a 
rare instance may sometimes occur, in which the length of 
time a comet continues visible, and a great number of obser- 
vations extending over a considerable portion of the orbit, 
may mark so great deviations from a parabolic motion as to 
lead to a tolerably exact estimation of the elliptic elements. 
For a single appearance we must be content with supposing 
the orbit to be a parabola ; a supposition which, if it be not 
rigorously true, serves important purposes in astronomy: it 
proves that the comets move round the sun by the same laws 
as the planets ; and it enables us to discover the identity of a. 
comet with one already observed, when we find that they agree 
in having the same parabolic elements. 

Three geocentric observations of the longitude and latitude 
of a comet are sufficient for determining the parabola which it 
describes. The problem is one of great difficulty. The appa- 
rent motion of a comet is the combined effect of its own motion 
and of that of the earth ; it is therefore extremely irregular 
and intricate ; and on this account it is difficult to deduce the 
heliocentric positions from observations made on the earth's 
surface. We can observe the planets at all times and in all 
situations ; and with regard to them we can thus select those 
positions where the heliocentric places are found immediately 
from observation, without any perplexed calculations; but we 
are deprived of this expedient in the case of the comets, 

Rs 
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which continue visible for a short time only, and in a small 
part of their orbit. 

In order to evade the difficulties attending a direct consi- 
deration of the problem, and to obtain an approximate solution 
at least, Sir Isaac Newton proposed to take a small portion 
of the orbit for a straight line described with a uniform motion. 
On this supposition the projections of the comet on the plane 
of the ecliptic will lie in one straight line as well as the real 
places in the heavens; and the -.several parts of both lines will 
have the same proportions as the intervals of time between 
the observations : so that in order to find the projection of a 
comet's trajectory on the plane of the ecliptic, we have only to 
draw a straight line which shall cut the several straight lines 
whose positions are determined by the observed longitudes 
in such a manner that the intercepted segments shall have 
given proportions. With three observations only this problem 
is indeterminate, or admits of innumerable solutions: for, let 
AP, BE, CQ, (fig. i. PI. IV.) drawn in the plane of the ecliptic, 
represent the directions in which the comet is seen from the 
earth at the three observations; then, having assumed any 
point at pleasure, as E, in one of those lines,* we can draw, 
through that point, a straight line PQ to terminate in the other 
two lines AP and CQ, so that PE and EQ shall have the same 
proportion as the intervals of time between the observations. 
When four observations are employed, the problem, gene- 
rally speaking, is determinate, and is easily solved.-f 

* Prin. Math. Lib. 3, Lem. 7. 

f Arith. Univ. Prob. 52. Prin. Math. Lib. 1, Lem. 27 Cor. Tho. Simpson's 
Elem. Geom. 3d edit. p. 353. Da v. Greg. Astr. Phys. et Geo. Elem. Lib. 5, Prop. 1 1. 
fidin. Tran. Vol. III. 
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If we reflect on the suppositions made by Newton for sim- 
plifying the method of finding a comet's trajectory, it must be 
admitted that they seem judicious, and that they are founded 
in fact. In the usual circumstances in which terrestrial obser- 
vations are made, there is no doubt that the real path of a 
comet, for a short time at least, differs little from a straight 
line, and that the motion is nearly uniform. Yet the problem 
grounded on these suppositions, although it has been often 
applied, has never, in any instance, led to a satisfactory deter- 
mination of a comet's orbit. The reason of this was first 
noticed by Boscovich, who shewed that in the actual state of 
the data, owing to the earth and the comet being both in mo- 
tion, the problem is as indeterminate when four observations 
are employed, as we have already remarked that it is in the 
case of three observations only. 

In order to explain what has just been said, it must be ob- 
served that, for a short time, the earth as well as a comet, may 
be supposed to move in a straight line with a uniform velocity. 
If we draw a tangent to that point of the earth's orbit which 
the planet would occupy at the middle instant between the 
two extreme observations, the supposition of a uniform motion 
in the tangent will not much displace the earth from its real 
positions, and will produce but little alteration in the apparent 
places of a comet. Suppose then that the earth moves with a 
uniform motion in the straight line AD, (fig. 2. PI. IV.) so that 
the parts AB, BC, CD, are proportional to the intervals of time 
between four observations of a comet : also let the four observed 
places of the comet in the ecliptic be represented by the points 
P, E, F and Q lying in one straight line of which the parts PE, 
EF, FQ, are proportional to the same intervals of time, and 
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consequently to the segments AB, EC, CD : then the lines AP, 
BE, CF, DQ, or the directions in which the comet is seen from 
the earth, will be given by position on account of the observed 
longitudes. Draw DM parallel to AP, and PM parallel to AD, 
and join QM : now if we assume any point whatever in the 
line AP, as/>; and draw pm parallel to PM, and mq parallel 
to MQ; then, the line pq being drawn, it will be cut by the 
lines BE and CF into segments having the same proportions 
as the segments of the lines AD and PQ. To demonstrate 
this : draw BH and CG parallel to AP and DM, and let those 
lines meet PM in H, G, and pm in h, g : join HE, GF ; and 
because the segments of PM are equal to the segments of AD, 
they will be proportional to the segments of PQ (hyp.); 
therefore (E. 2. VI.) HE and GF will be parallel to MQ and 
mq: drawee and gf parallel to HE and GF; and consequently 
(E. 2. VI.) pq will be cut in e and/ in the same proportion as 
AD is cut in B and C, and PQ in H and G: and I say that BE 
will cut pq in e, and CF will cut it inf. For on account of the 
parallel lines, and because PM =pm, PG — pg, PH = ph, 
therefore 

QM : HE : : qm : he 

QM-.GF ::qm:gf 
and, by alternation, 

QM : qm: : HE : he 

QM : qm : : GF : gf 
but, on account of the similar triangles QDM, qdm, 

QM : qm : : MD : mD 
therefore, because MD = GC = BH, and mD = Cg = Bh, 

BH : Bh : : HE : he 

CG:C g :>GF:gf 
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consequently BE will cut pq in e, and CF will cut it inf. Thus 
it is proved that innumerable lines may be drawn that will all 
be cut by the four lines AP, BE, CF, DCj into segments hav- 
ing the same proportions as the intervals of time between four 
observations of a comet. The problem is therefore not more 
determinate with four, or even a greater number of, observa- 
tions, than it is with three : and it is totally unfit for finding 
the lengths of the lines AP, BE, CF, DQ, or the comet's dis- 
tances from the earth. The shorter we make the intervals 
of time between the observations, and the more exactly the 
suppositions made by Sir Isaac Newton are fulfilled, the less 
fitted is the problem to answer the purposes of astronomy, 
because it approaches more nearly to the conditions which 
make it indeterminate.* 

Besides the method of which we have hitherto been speak- 
ing, Sir Isaac Newton likewise gave another very able solu- 
tion of the problem,"}- which however is too tedious and 
laborious in practice to fulfil the wishes of astronomers, and 
it is now disused. Dr. Halley, who studied the astronomy 
of comets with much diligence, has no where explained the 
method he followed in determining their orbits. Since his time 
this problem has not been much discussed in England ; but it 
is the subject of a great number of researches scattered in the 
academical collections of the learned societies on the conti- 
nent. All the resources of the modern mathematics have 
been directed to overcome the peculiar difficulties of this in- 
tricate investigation. Yet it must be confessed that practical 
astronomy has not reaped so much benefit as might justly 

* Vide a Memoir by Professor Playfair in the 3d Vol. of the Edin. Trans. 
f Prin. Math. Lib. 3, Prop. 41. 
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have been expected from the labours of so many eminent 
men. In most of the solutions that have been proposed we 
find certain coefficients which are not only extremely small 
when reduced into numbers, but they likewise suffer great 
alterations in their relative magnitudes when small changes, 
within the limits of the probable errors, are made in the quan- 
tities furnished by observation. In such circumstances the 
results obtained are often wide of the truth, and in no case can 
entire dependence be had on their accuracy. The fact is that the 
solutions here spoken of, are uncertain and useless in practice 
for the very same reason as the one already mentioned, which 
supposes the orbit to be rectilineal, and proposes to deduce it 
from four observations Nor is this to be wondered at. We 
have already shewn that it is not any inaccuracy in the sup- 
positions made by Sir Isaac Newton, which renders his first 
method insufficient and of no use in practical astronomy; nor 
is this owing to any defectiveness in the geometrical construc- 
tion ; it arises from his having overlooked that connection 
between the motion of the earth and the motion of the comet, 
which makes the data approach so near the indeterminate case 
of the problem that the conclusion becomes quite uncertain. 
We fall upon the indeterminate case when we suppose both 
the earth and the comet to move in straight lines with uniform 
velocities; and the very same hypothesis will be found to 
make the small coefficients alluded to, accurately evanescent.* 
The inconvenience is inherent in the quantities obtained by 

* The small coefficients here mentioned, are exactly evanescent when all the three 
observed places of a comet are in one great circle of the heavens : but it is easy to 
prove that if two celestial bodies be supposed to move uniformly in straight lines, one 
of them will be seen from the other to describe constantly the same great circle. 
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observation, and it cannot be entirely obviated merely by 
pursuing new methods of solution, nor by pushing calculation 
to greater degrees of exactness. The problem cannot be 
justly and usefully solved, unless such coefficients are ex- 
cluded, and the result is obtained by means of quantities on 
which the errors of observation have no more influence than 
they ought to have. 

There are however some solutions of this problem to which 
the preceding observations must not be applied. Of this kind is 
the method of M. Boscovich ; that of the celebrated Laplace ;* 
and those which Legendre has more lately published : -f all 
of which have been found useful in practical astronomy. The 
method of Boscovich owes its utility as an approximation to 
the circumstance of introducing the velocity in the orbit as a 
principal condition: for that velocity depending upon the pro- 
portion of the distances of the earth and the comet from the 
sun, limits the other conditions, and places the orbit in its pro- 
per situation. The same thing may be said of the methods of 
Laplace and Legendre : and, in general, we may affirm that 
no solution of this problem can be free from the imperfections 
we have pointed out, in which the velocity in the orbit, or 
some equivalent property, does not enter as a principal con- 
dition. 

In order to place what has been said in a clearer light, it is 
to be observed that three complete observations of a celestial 
body are sufficient for determining the species, the magnitude, 
and the position of the curve in which it moves round the sun. 
On this account there is a superfluity of conditions when we 

* Mec. Celeste, Tom. i. Liv. z, chap. 4. 

f Nouvelles Methodes pour la Determination de» Orbites des Cometcs, i8c6. 
MDCCCXIV. S 
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suppose the orbit to be a parabola : because, in this ease, the 
velocity in the orbit furnishes an equation without introducing 
any new unknown quantity. Thus it happens, that in the 
problem of the comets there is one equation more than there 
are quantities sought: and by combining those equations in 
different ways, various solutions of the problem may be ob- 
tained. But it ought likewise to be observed, that if we set 
aside the equation derived from the nature of the orbit, the 
remaining ones, in the actual state of the data, will nearly 
coincide with what would result from the hypothesis of a 
uniform motion in a straight line: and although, theoreti- 
cally speaking, we can solve the problem by means of those 
equations, yet we shall thus infallibly introduce coefficients 
that are small and ill defined, and unfit for any practical 
purpose. It is therefore necessary to include the velocity in 
the orbit, or some equivalent property, if we wish to obtain a 
solution useful in practical astronomy : and even when this 
mode of solution is adopted, it is still necessary to examine 
with care the quantities introduced by combining the other, 
conditions, in order to exclude the faulty coefficients we have 
been speaking of. 

It is of the greater consequence to discuss the peculiarities 
of this problem, because the observations of comets are sus- 
ceptible of little accuracy even with the best instruments and 
the greatest care, on account of the haze, or coma, with which 
those bodies are generally surrounded. Every solution of a 
physical problem which is deduced from quantities that are of 
the same order as the unavoidable errors in the data furnished 
by observation, can plainly be of no practical utility ; and it is 
in this predicament that those coefficients stand, which would 
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vanish exactly on the supposition that the earth and a- comet 
moved, for a short time, with uniform velocities in straight 
lines. 

On the whole it appears from What has been said, that the 
same circumstance to which is owing the failure of Sir Isaac 
Newton's first method of finding a comet's trajectory, has 
produced like effects in many of the later solutions of this 
problem. This is indeed no more than what was to be ex- 
pected : and it must excite surprise that more attention was 
not paid to the peculiarities of this investigation, after Bosco- 
vich had fully developed the solution of Newton, and shown 
the cause of its want of success. 

The English astronomer will find the methods of Boscovich 
and Laplace, properly illustrated by examples of all the 
computations, in a work published in 1793, by Sir Henry 
Englefield : in which the author has judiciously selected all 
that is practically useful from the numerous writings on this 
subject at that time before the public. 

The method of Laplace would be rigorous, were it possible 
to find exactly the numerical quantities that enter into his for* 
mulas. These are the first and second differential coefficients 
of the longitude and latitude considered as functions of the 
time; the values of which we can do nothing more than deter- 
mine nearly by interpolating the observed places of a comet. 
With three observations only, the quantities thus found can- 
not be expected to have much accuracy, more especially if 
the motion in longitude, or latitude, be quick and irregular, as 
often happens; and when a greater number of observations is 
employed, there arises a new cause of inaccuracy in the arith- 
metical operations necessary for interpolating, which augment 

S 2 
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the errors of observation, always considerable of themselves 
in the case of the comets. These remarks, which are cer- 
tainly just, are made by Legendre in his memoir on this pro*- 
blem, published in 1806. That geometer therefore thinks it 
preferable to deduce the orbit immediately from the observed 
places of a Comet: and we agree with his opinion, that the 
observations will be better represented this way, than by em- 
ploying the method of Laplace. With regard to the solutions 
of the problem which Legendre himself has given, it will be 
admitted that they approximate sufficiently near to the ele- 
ments of the orbit, to answer all the purposes of practical as- 
tronomy ; but his formulas are complicated ; and the number 
of equations which it is necessary to form and to solve, seems 
to render his methods ill adapted for general use. It is the 
object of this paper to give a new solution of this problem, 
which, while it does not yield to any of the known methods 
in accuracy of result, I judge, will be found as commodious in 
practice as the nature of such a calculation can well admit. 

1. Let the coordinates that determine the position of a 
comet with regard to the ecliptic, and a straight line drawn 
through the first point of Aries, (which point is supposed to 
be immoveable) be represented by x,y, z; of which the last 
is perpendicular to the plane of the ecliptic, and the other two 
have their origin in the sun's centre : and likewise let the 
comet's distance from the sun, or the radius vector of the or- 
bit, be represented by r. Farther, supposing the mean dis- 
tance of the earth's orbit to be unit, put m for the length of 
the circular arc of the mean motion described in a second of 

time; then —, the versed-sine of that arc; will bef the space 
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through which the earth descends towards the sun in a second 
of time; and m* is the velocity produced by the solar force; 
therefore the direct velocity towards the sun produced by the 

same force acting on the comet at the distance r, is -j ; and, 

this being decomposed into the equivalent velocities parallel 

to the coordinates, these velocities will be — , ^-£, ~ re - 

T T T 

spectively, which are the actual diminutions of the velocities 
with which the coordinates increase ; but if dr represent the 
fluxion of the time, the diminutions of the velocities with which 

the coordinates increase, will be represented by — -p-, — ~> 

— ~; therefore we shall have these equations, viz. 

m*x ddx 4 m*y ___ ddy # m*x _ ddx 

">"' ~" dF* ri 17*' ~t* IF' 

Now if we put t to denote the mean motion of the earth, 
proportional to the time, and estimated by the arc of the circle 
whose radius is the mean distance, or unit ; then dr* will de- 
note the same thing that is signified by mVr 2 in the foregoing 
equations, which will thus become, 

ddx , x 

7? + H = ° 

# + * — (X) 

ddx. , x ^ 

IF "T IF — °' 
Multiply these equations by zdx, sdy, zdz, respectively; then, 
because zxdx + zydy + szdz = zrdr, we shall get by adding 
them 

2dxddx+ 2dyddy + zdzddx , zdr 

j? r r * — ; 

and, by integrating, 

g±j£t£ - i « const. 
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To determine the constant quantity in this equation, it is to 
be observed that the first term is the square of the velocity in 
the orbit : if then we suppose the orbit to be an ellipse of 
which/ is the mean distance ; f. t the eccentricity ; and u and t/ 
the velocities at the perihelion and aphelion distances; the 
foregoing equation will become at those two points of the 
orbit, 

w a — ttt—x = const. 

u ' ~7oTD = co " st - 
and if we multiply these by (1 — if and (l + 0* respectively, 
and then subtract them, we shall get 

u 2 (i — s y — v n .(i + 2 + j= — ¥ x const. 

but u ./( l— <•) — xt'.f ( l 4- < ) = o ; for these quantities are 
respectively the doubles of the sectors described at the peri- 
helion and aphelion in the time denoted by unit : hence j 

= const. : the foregoing equation will therefore become, 
dx*+dy*+dz* 2 , i 

It* ~~ r T / — y - 
Again let the equations ( i), multiplied by x,y, z respectively, 
be added to the last equation ; then observing that xddx + 

yddy -j- zddz -f- dx* -f- dy* -\- dz 2, — — ~— , we shall get 

i — r * i , i 

a" * dr» *~" r "•" 7" — °' 

ft r* 

Let r° denote the value of r, and 2/> that of ■—, at the epoch 

from which the time is reckoned ; and further let the value 
of r* be assumed in a series with indeterminate coefficients, as 
follows, viz. 

r- = r" 2 + 2/>r + At 2 + Bt 3 &c. : 
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then — = ~ &c. ; and if these values be substituted in 

the foregoing equation the coefficients A, B, &c. will be readily 
determined ; which being done, we shall get 

It is next necessary to integrate the equations ( 1 ). I begin 
with seeking an incomplete integral of this form, viz. 

s = 1 -f Ct* + Dr 3 &c. : 
and, having substituted in the first of the equations (1), s for x; 
and, for -^5, its value taken from the preceding series for r a ; 
I get, 

by means of which equation, the coefficients in the series for 
s y are easily found: then 

*=!-*• 5 + -f.£&C (*)• 

I now seek another incomplete integral of this form, viz. 
0- = r + E . t 3 + F . t 4 &c. : 
and, having substituted <r for x, I get, 

£ + £•{«-£**}=>.: 

whence the coefficients of the series for a- will be readily 
found; then 

r.= T-*..£ + .£.£&c (3). 

The complete integrals of the equations (1) can now be 
obtained, and they are 

c- , dx- 
X — SX + v -% 

o 1 dy 

y — sy + <r-£ 
% = s% + <r ^ : 
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where x° and —-, which denote what x and ■£ become at the 

epoch, are the arbitrary quantities necessaryto complete the 
fluent ; and the same things are to be understood in the ex- 
pressions of y and z. 

2. Of three observations which are necessary to determine 
a comet's orbit, let the epoch, or the instant from which the 
time is reckoned, be made to fall on the middle one ; and let 
— t denote the interval between the first and second obser- 
vations, which precedes the epoch ; and -f- t', the interval 
between the second and third observations, which follows the 
epoch ; and further let x, y, % and r express the coordinates 
and radius vector of the comet at the first observation ; x°,y°, 
%*, and r°, the same quantities at the second observation ; and 
x\ y', z', and r', those at the third : finally, let the values of 
the series denoted by s and a- (Equat. 2 and 3) be taken for 
the intervals — t and -f- r ; that is, let 
_„ 1 T * P* 3 

S - — * 3" • 7°» ~zr»s 

t 3 1 P** 



IX 

j-oS I 2J" oS 

6 " *°* ' 4r«» ' 






then, on account of the equations (4), we shall get, 



</*° 

iT sn SX *— • (T -j— 
Z =z SZ — a- -^ 

x = *'■* + "• "37 
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and if these equations be combined two and two so as to ex- 
terminate ■—, ~^, ~,we shall get 

cr'x -f- vx 1 = ( °"' 5 + °~ s ' ) x * 

Sy + *?=*{* + **)?. (5) 

<t% -f" <r%' = (<r's + its') X* : 
Finally let <r" be written for (<r's -J- «') ; then, by substituting 
the series that <r, </, 5, s' stand for, we shall get 

§. Employing the same marks as in the case of the coor- 
dinates, to distinguish the known quantities of each observation, 
I shall write, 

c, c°, (f ; for the three geocentric longitudes of a comet : 

x, x°, x' ; for the three geocentric latitudes : 

p, p°, p ' ; for the comet's distances from the earth : 

e, e°, e'\ for the earth's longitudes : 

R, R°, R'; for the earth's distances from the sun : 
by means of which quantities, the coordinates of the comet at 
the first observation will be thus expressed, viz. 
x = R cos. e + p cos. x cos. c 
y = R sin. e -\- p cos. x sin. c 
z = p sin. x : 
and in like manner the expressions of the coordinates at the 
second and third observations will be formed by placing the 
characteristic marks on the several letters. By substituting 
the values thus obtained in the equations (5), we shall get 
these other equations, viz. 

mdcccxiv. T 
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<r'p cos. X cos. c -{- <rp' cos. x' cos. c' — ^"p" cos. x° cos. c" = 
— ar'R cos. e — o-R' cos. e' + ^"K 3 cos. e° 

<r"p cos. x sin. £ -J- °~p' coS - x> s " 1 * e ' ~H °"V cos - x ° sm - c ° == (7) 

— ^'R sin. * — o-R' sin. e' -f <r"R° sin. <?° 
<r'p sin. x + °"p' sin. x' — <r"p° sin. x° = o. 

These equations are universally true of all orbits ; and, in 
order to solve the general problem of finding a planet's orbit, 
they require only to be properly discussed with the view of 
obtaining the values of p, p", p' by means of formulas sufficiently 
exact and commodious in practice: but, at present, I confine 
my attention to the case of the comets moving in parabolic 
orbits. 

4, We may apply the equations (5), which are general for 
all orbits, to the orbit of the earth : let jtt, /*', p" represent 
what <r, <r', a-" become in the case of that orbit; that is (neg- 
lecting quantities of the fourth and higher orders), let (Equat. 

s, 3> 6 ) 
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then, Equat. (5), 

p'R cos. e 4- fiR' cos. <?' — |a'R° cos. e° = 
pR sin. e -f- [*R sin. e' — p'R° sin. e" = 0, 
On account of these formulas, we have 
tr'R cos. e 4- <rR' cos. e' — a"R a cos. e" = 
(<r'~ {//) R cos. <? + («■ — p) R cos. <?'_(</'—/') R°cos. e°: 
and, by substituting the series that <r, <r', <r"; p, /*', /a"; stand 
For, in the right hand side of this equation, and neglecting 
quantities of the fourth and higher orders, we shall get 
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<r f R cos. e -f <rR' cos. e' — <r"R° cos. e" = 
T'&--£*) ■ {-t 3 Rcos.,-t 3 R'cos.,'+(t4.t') 3 R cos., } 
or, which is the same thing, 
<r'Rcos.e-f<rR' cos.e'— <r"R° cos.e°— !^±^ . (_L„_J_).R° C0S , (? » 

. — T * (£» — 55) • ( R cos - e ~ R ' C0Si e °) 

- T • (i ~ Es) ' ( R ' cos « e '~ R ° cos - O •■ 
but (R cos. e — R° cos. <?°) and (R' cos. e' — R° cos. e°) are 
quantities of the first order ; consequently the two terras that 
contain them, will be of the fourth order, and may therefore 
be rejected : thus we get 

<r'R cos. e + a-R' cos. e'— <r"R° cos. e° = ££±2 . |i- — J- f J . 
R° cos. <?°. 
And, in the very same manner it may be shewn that 

(r'R sin. e + o-R' sin. <?' — <r"R° sin. e° =» tllil±l2 . |_L _ J_ j . 
R° sin. e°. 

Let the values just obtained be substituted in the two first 
of the equations (7), and they will become 
«fp cos. x cos. c -J- 07/ cos. x' cos. c' — <r"p° cos. x° cos. c° s= 

<r'p cos. x sin. c -j- cp' cos. x' sin. c' — tr"p° cos. x° sin. c° =s 

- ^M^« " ^) • R ° sin - *° : 

Further let n denote an angle, to be afterwards determined ; 

then, by combining the two last equations we shall readily 

obtain, 

ff-'p cos. x cos. (c — n) -f- 07/ cos. x' cos. (c f — n) — <r"p° cos. x° 

x cos. (c°- n) = - n^ . (i - £) . R° cos. (^) 

Ta 
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o-'p cos. x sin. (c — ») + rf cos - x ' s i n - C^— n ) •— ^V cos. x e 

x sin. (<•-*) = - ^±$ . (£ - i) • R ° *«• (*°- »)• 
Conceive a great circle to be drawn on the surface of the 
sphere through the geocentric places of the comet at the two 
extreme observations; and let #, in the foregoing equations, 
denote the longitude of one of the intersections of that great 
circle with the ecliptic: also let h and h' denote the arcs of the 
great circle intercepted between the same intersection and the 
places of the comet abovementioned ; then these arcs will be 
the hypothenuses of two right-angled triangles of which the 
sides, perpendicular to the ecliptic, are the arcs x and x'; and 
the other sides, in the ecliptic, are the arcs (c— n) and (cf—n) : 
further, let i denote the inclination of the same great circle to 
the ecliptic; then will i be an angle common to the two right- 
angled triangles abovementioned, opposite to the sides x and 
x'; and, by the properties demonstrated in spherical trigono- 
metry, we shall get these formulas, viz. 

cos. x cos. (c — n) = cos. h 
cos. x' cos. (c' — n) = cos. k' 
cos. x sin. (c — n) = cos. i sin. h 
cos. x' sin. (c' — n ) = cos. / sin. h' 

sin. x = sin. i sin. h 

sin. x'= sin. i sin. h'. 

Now let these values be substituted in the two foregoing 
equations, and likewise in the last of the equations (7), and 
we shall get, 
(r'p cos. h -f- <rf>* cos. h' — <r"f cos. X° cos. (c° — n) = 

-^•(i-^)-* «■ ('-») (8) 
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<r'p cos. i sin./* -f °"p' cos « * sm - h'—<r"p° cos. x° sin. (e — n) = 

-*£*2. ■(£-£). IP *>.('-») (8) 

o-'p sin. z sin. /i + ff p' si n » * sin. '^'- — 0""p° s ' n - x ° — 0: 

and finally if we subtract the two last of these equations after 

having divided them by cos. i and sin. i, we shall get 

// o I cos.?.°sin . (<: — n) sin A* ^ (t+t')tt' fi _i_\ R*sin (<°— n) ^ . 

° 1 cos. i sin. i J 2 \r° 3 R o3 j " cos. i Wl* 

5. In the last equation, all the terms being of the third 
order, we may consider <r" as equal to (t + *"') : therefore, by 
division, we shall get 

o cos. x» f . , v tan. K* 1 tt' / I I \ „<, 

P *sTnT(7C7r l Sin -( C - M )--EnTTS == T-lr--R^!- R i 

or, by introducing a new letter, we shall have these two for- 
mulas, equivalent to the last one, viz. 

c- COS. >.* f • / o ^ tan< x ° 1 

£_. ___. — ^ . ] sin. (c— n ) — - — r t 

* sin. (e°-») C V I tan. iJ fir>\ 

5. o tt'. / I I ] po 

If we suppose all the three geocentric places of the comet 
to be situated in the same great circle of the heavens ; then 
sin. (c°— n) . tan. i = tan. x°: in this case therefore £= o, 

which requires that ^ — ~ = o, and r° == R°. Thus we 

learn that a comet and the earth are equally distant from the 
sun, when the comet's apparent motion continues for a short 
time to be performed in one great circle of the heavens. Again 

if £ be positive, then (— — ~) will be positive too 5 and r° will 

be less than R°: but if £ be negative, then f ^ — g-J will like- 
wise be negative, and r° will be greater than R°. But there are 
two cases when £ will be indeterminate, and the preceding rules 
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will fail. This will happen when the comet moves exactly in 



the plane of the ecliptic, in which case — -'— ~ is indeterminate. 



tan. I 



The same thing will likewise happen when the great circle 
which passes through the two extreme places of the comet 
cuts the ecliptic in the points occupied by the earth and the 
sun at the middle observation : for in this case, sin. (e° — n) 
= o ; now since £ cannot become infinitely great, it is necessary 

that {sin. (c° — n) — 7^-7} vanish together with sin. (e°~- n), 

which will make £" indeterminate. 

The formula for £", which by its sign alone enables us to 
judge of the magnitude of r°, contains the angles n and i : and 
it is necessary to investigate rules for computing those quan- 
tities. By attending to what n and i were made to signify it 
will readily appear that 

tan. a. tan. ti 

tan. I =• -r— r = -: — ~, i : 

sin. (c— n) sin. (c — n) 

let /=lnf75 then 

/sin. {c — n) = sin. {c' — n): 
but 

/ . fc'+* % c' — c lc'+c , . c'—c 

sm. {c— n) = sin. l-~ n) cos. — cos. i-j n) sin. - — - 

'it \ • lc'+c , c — c . Ic' + c » . c'—c 

sm. [c —n)sx sm. \-~ «)cos. -j- -j-cos. l-j n)sm. i —~ i 

therefore by substitution and division we shall get 

lc'+c n f+l , c'—c 

tan. -J n) =j^~ t . tan. — : 

and, if/ = £j^ = tan. (45 ± <p) ; then 

tan. (*±£ - n) = ± ^S = + tan. £=S cot. <p. 
To the value of tan. (~^ — n), found by this rule, there 
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will correspond two arcs, differing from one another by 180 , 
which will determine both the points of intersection of the 
great circle and the ecliptic ; but it will be sufficient to take 
that one which immediately precedes the comet in longitude. 
The arc n being thus found, we have these formulas for 
finding i, viz. 

tan. x tan. x' 

tan. t = 



sin. (c— n) sin.(c— «)' 

and the double value will serve to prove the accuracy of the 
calculation. 

But we may determine whether r° is greater or less than 
R°, without any calculation. Let h° denote the arc of a great 
circle drawn from the intersection whose longitude is n, to the 
geocentric place of the comet at the middle observation : then 
h° will be the hypothenuse of a right-angled triangle, having 
the arc (t° — n) of the ecliptic for one side, and x°, perpendi- 
cular to the ecliptic, for the other side: and if we put i° to 
denote the angle of the triangle opposite to x°, we shall have 
cos. x° sin. (t° — n) = cos. i° sin. h* 
sin. x° = sin. i° sin. h°. 
Substitute these values in the equation (9), and we shall get, 
by division, 

o sin. (i—i a ) ■ sin &" tt' /j_ __ _J_\ tio 

P ' sin. i sin. («°-w) — 2 ' [r>* Ro>j • R ' 

Conceive two arcs to be drawn perpendicular to the great 
circle which passes through the geocentric places of the comet 
at the two extreme observations ; one, denoted by tt, drawn 
from the extremity of If, that is, from the geocentric place of 
the comet at the middle observation ; and the other denoted 
by n, drawn from the place in the ecliptic occupied by the 
earth at the same time: then, observing that it will be one 
side of a right-angled spherical triangle, of which h° is the 
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hypothenuse, and (i — i ) the angle opposite to n ; it is plain 
that we shall have 

sin. it = sin. (i—i°) sin. h° 

sin. n = sin. i sin. {e°— - n). 
Therefore, by substitution, the last equation will become 

o sin.gr t£ ll_ _i_\ Ro 

P x sin. n ~~ 2 .' \r° 3 RoJJ • ^ * 

This equation is equivalent to the second of the equations 
(io): and we infer from it that r° will be less, or greater, than 

R°, according as j|~ is positive or negative ; that is, accord- 
ing as the arcs it and n are on the same side, or on different 
sides, of the great circle to which they are both perpendicular. 
Hence if we mark, on the surface of a celestial globe, three 
geocentric places of a comet, and likewise the places in the 
ecliptic occupied by the earth at the middle observation ; then 
the comet's distance from the sun will be greater than the 
earth's distance, when the great circle drawn through the two 
extreme places of the comet passes between the earth's place, 
and the remaining place of the comet ; but when these two 
places are both on the same side of the great circle, the comet's 
distance from the sun will be less than the earth's distance. 
When all the three places of the comet are in one great circle 

of the heavens, then <n = o, and s ^~ = o : and in this case, 

the comet and the earth will be equally distant from the sun. 
From these rules we must however except the cases in which 

sin. n = o : for, since ^~ cannot become infinite, it is neces- 
sary that the numerator vanish together with the denominator, 
so that 2575 will become indeterminate. This will happen, as 
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before observed, when the comet's orbit coincides with the 
ecliptic, in which case sin. i = 0; and likewise when the great 
circle drawn through the two extreme geocentric places of the 
comet cuts the ecliptic in the places of the earth and the sun 
at the middle observation, in which case sin. (e° — n) = 0. 

The inferences here drawn from the preceding analysis 
coincide with the rules first given by M. Lambert, of Berlin, 
forjudging of a comet's distance from the sun by the inflection 
of its apparent path. 

6, The second of the equations (10) contains only one un- 
known quantity, namely p° : and hence it may be thought that 
we have already, by means of that equation, obtained a solu- 
tion of the problem, which is both simple and elegant. And 
this would undoubtedly be the case, were it not that the co- 
efficient £ is always extremely small and greatly affected with 
the errors of observation. It depends entirely on the deviation 
of the comet's apparent path from a great circle of the heavens; 
and this deviation is often so little, that small changes in the 
observed places of the comet, by no means inconsistent with 
the errors of observation, will make £ evanescent, or even 
take a different sign from what it had before. If we suppose 
the motions both of the earth and the comet to be rectilineal 
and uniform, which is never far from the truth in the short 
intervals that must intervene between the observations selected 
for finding a comet's orbit; then the apparent motion of the 
comet would be accurately in a great circle of the heavens, 
and the mode of solution here alluded to could not be applied 
at all. 

The equation here spoken of may no doubt be usefully ap- 
plied in favourable circumstances, particularly when the comet 

mdcccxiv. U 
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is at a considerable distance from the sun, and has a slow mo- 
tion, allowing an interval of 20, or even a greater number of 
days, between the extreme observations ; in which time the 
apparent deviation from a great circle of the heavens will be 
considerable enough to define the magnitude of £ with tole- 
rable certainty. But in many circumstances the same mode 
of solution cannot be applied with advantage or success : as 
when a comet has a quick motion, which will allow only a very 
short interval between the observations ; or when a comet 
and the earth are both at the same distance from the sun, or 
nearly so, in which case £is evanescent, or very small. 

It appears then, that, generally speaking, the problem of 
the comets must, in reality, be considered as indeterminate, if 
we set aside the condition derived from the nature of the pa- 
rabola. In this respect there is an essential distinction between 
the investigation of the orbit of a planet and that of a comet ; 
a distinction which, when due regard is had to practical utility, 
renders it necessary to separate the particular case from the 
general problem. When the question is to determine a pla- 
net's orbit by means of three observations made at short in- 
tervals of time, without any assistance from the supposition 
of a circular orbit, or other hypothesis, the equations which 
the data furnish are barely sufficient for finding three radii 
vectores of the orbit and the two angles contained between 
them; and there are no supernumerary conditions to compen- 
sate what may be imperfect in the observations. In whatever 
manner we proceed, the solution will lead to an equation, such 
as that we have been speaking of; and the success of the inves- 
tigation will therefore depend upon our being able to ascertain 
with some degree of certainty the magnitude of the coefficients 
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which that equation contains. If there were the same uncer- 
tainty in this case that both theory and experience agree in 
shewing there is in the case of the comets, the problem, as 
applied to the planets, would be a mere theoretical speculation 
without any practical utility. It may therefore be asked, what 
are the precise circumstances which make the investigation 
successful in the one case, while in the other we must seek 
aid from the nature of the orbit in order to obtain a useful 
result? In the first place, the observations of the planets are 
susceptible of a much greater degree of accuracy than those 
of the comets ; which confines the uncertainty arising from 
the errors of observation within much narrower limits in the 
one case than in the other. In the second place, the distance 
of a planet from the sun is never equal, or very nearly equal, 
to the earth's distance from the same body ; and this occasions 
a greater curvature of the apparent path than often takes place 
with regard to the comets. We may likewise add that the 
motion of a planet, in an orbit of moderate eccentricity, is 
slower than the motion of a comet at a like distance from the 
sun ; which allows a longer interval between the observations 
of a planet, and thereby contributes to heighten the effect 
arising from the inflection of the apparent path. It is to these 
causes, as well as to the excellence of his methods of investi- 
gation, that we must ascribe the great exactness with which 
the orbits of the new planets discovered in the course of the 
present century, have been determined by M. Gauss, in a work 
of great merit, on the Theory of the Motion of the Celestial Bodies, 
published in 1809, which cannot but add much to his reputa- 
tion, already very high on account of former scientific disco- 
veries. 

Us 
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7. If we put 

» 1 S • / o \ tan - x ° \ 
$ = - — r -r— r . « sm. (c° — «) — t \ ; 

tan, («*—») (.•■..* ' tan. z J ' 

then the equation (9 ) will become by substitution 

«r"cos. x«; J. ,•« ££±2 . (i - pi) • R° cos. (e%-»): 

and this value being substituted in the first of the equations 
( 8 ) , we shall get 

<r'p cos. h -}- <rp f cos. h' =sz-<r"p° cos. X° . {cos. (c°— n) — Si. 
But, by putting cos. x° tan, x° for sin. x°, and then dividing by 
sin. i, the third of the same set of equations will become 

o-'p sin. h -{- <rp' sin. W = <r"p° . cos. x° . — v , 
Now if we make 



1 tan. W 3 

tan. u = t — : x 



sin. i cos. (c°— «) — £ ' 

then, by subtracting the two equations just found, after having 
multiplied them by sin. a and cos. a respectively ; we shall get 

r'p sin. (w — h) — <rp' sin. (h! — «)=o 
whence 

f sin. («— b) a-'' 

P P X sin (b — u) ' T* 

The value of — , obtained by substituting the series that </ 
and w stand for, will be as follows, viz. 

£_ z_ S fjzzll 

a- T * l 1 6ro3 i' 

In this formula -~ 6 I- is evanescent when t = t'; and in all 

cases it is a very small quantity; because the intervals of time 
between the observations made use of for finding a comet's 
orbit, should in no instance be extremely unequal. If then 

we suppose — • = -~, we shall get these two formulas, viz. 
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n — sin; ( ".-"ft) t ' 

p'=/3. P ; ^ ' 

which are very exact, when the intervals between the obser- 
vations are not too great. 

The angle « which enters into the preceding formulas, de- 
pends partly upon the value of the quantity $; and it thus 
becomes impossible to assign the value of that angle when $ 
cannot be determined. If we compare $ with the quantity 

formerly denoted by £(No. 5), we shall get £= £*x cos J * ^* ■ : 
hence it appears that £ will be indeterminate when £is so; and 
in such cases therefore the last formulas will fail. Now the 
cases in which £ becomes indeterminate have already been 
noticed (No. 5); they differ by a real distinction from the 
other cases of the problem, and require a separate discussion, 
if we wish to have clear and precise notions on the subject of 
this research : we shall therefore return to the examination 
of them in the sequel. 

8. On account of the formulas (11), the coordinates of the 
comet at the two extreme observations, will depend only upon 
one unknown quantity, namely p : and if we substitute the 
values of the coordinates in the following expressions, viz. r a 
= j* 4. y 4. z * ; r 1 * = x n + /» + z n ; V = xx' + yy' + %%' ; 
there will result these other expressions, which contain no 
unknown quantity but p, viz. 

r 2 — R 2 + 2R cos. x cos. [e — c) . p -f- p* 
r" =■ R" + 2R'/3 cos. x' cos. (*'— d ) . p + /3 2 p* 
V = RR' cos. (*'— «r) + { R/3 cos. x' cos. (e — <;'■) -f R' cos. * 

X cos. (t-' — c) I X p + /3 COS. y . p 2 

cos. ^ = cos. X cos. x' cos. (c' — v) 4- sin. x sin. x'. 
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We must now have recourse to the properties of the orbit 
to get such an equation between the functions, r*, r n and V as 
shall serve to determine the unknown quantity which they 
contain. For this purpose we may employ the expression 
which gives the time of describing an arc of a parabola, by 
means of the chord of that arc and the two radii vectores drawn 
to its extremities. This elegant property of the motion in a 
parabola, seems to have been first found out by Euler, as 
M. Gauss informs us ; but it is commonly attributed to M. 
Lambert of Berlin, who probably discovered it without know- 
ing what Euler had done, and who has extended it to all the 
conic sections.* It would be superfluous to give here the in- 
vestigation of a truth so well known : it will be sufficient to 
refer to the Mecanique Ctteste of Laplace,-^ or to the work of 
M. Gauss. Let b denote the chord drawn between the places 
of the comet at the two extreme observations : then observing 
that t -|- r' is the arc of the earth's mean motion correspond- 
ing to the time of describing the parabolic arc of which b is 
the chord ; we shall get 

T -f-T'=|. {(r + r'+ft) 1 — (r + r'— bf}: 
and, by expanding the radicals, 
r +•= v/F+F. JA _ ± . JL-. _ -L . _£_ &,.}: 

and, by squaring and omitting the sixth and higher powers 
of b, 

Let a 2 = 2r* + 2r" = (r + r*) m + (r — r'f = 

* Theor. Mot. Corp. Codes. Lib. I. Sect. 3. § \o6. 
f Prem* Part, Lib. II. Chap, 4. § 27, 



a first Approximation to the Orbit of a Comet. 151 

(r+n' + t&^f- the " 

r + ,= 8 .{ 1 _i.(^- i .(^r & c.}: 

and by substituting this value in the foregoing expression we 
get, nearly, 

consequently 

which is exact to quantities of the sixth and higher orders. 

Now b* = (x-x'Y + {y — y'f + {z-z'f = x 2 +/ + z* 
- 2 (xx' + yy' + zz') + x n + y ,% + z' 2 = r 2 + r n - 2 V : 
therefore in this last equation all the quantities concerned de- 
pend upon r 2 , r n , and V, which themselves contain only one 
unknown quantity, namely p : thus that equation will serve to 
determine p, or the comet's distance from the earth at the time 
of the first observation, on the finding of which the solution 
of the problem depends. 

9. If we now collect in one view all the formulas that have 
been investigated, we shall have the following method for 
computing the distance of a comet from the earth. 

1st. The symbols 1 and r' denoting the intervals in days of 
mean time, between the middle observation and the first and 
second observations respectively ; we begin with computing 
the arc of the earth's mean motion in the time (t-\-t')\ this 
arc we shall now denote by 0, and it will be computed by this 

formula, viz. S = *? itg /g "» 2ar being the circumference of a 

circle whose radius is unit, and 365.25638 the length of the 
sydereal year ; so tnat the logarithm of will be found by 
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adding the constant logarithm 8.2355814 to the logarithm of 
t + t'. 

We must next compute the angles n and /, by the formulas 
investigated in No. 5: n being the longitude of one of the 
points in which a great circle drawn through the two extreme 
places of the Comet cuts the ecliptic ; namely, the longitude 
of that intersection of the two circles which immediately pre- 
cedes the comet in longitude : and i being the inclination of 
the same great circle to the ecliptic. 

We must further make these computations, viz. 
cos. h = cos. x cos. (c — n) 
cos. h' = cos. x' cos. (c' — n) 

S=z -= — ? . 1 sin. (c° — n) — - — r-\ 

— tan. («•—«) L v ' tan. i i 

1 tan. x° 

tan. u = ~i — , x 



sin. i cos. (c* — n) — 8 
„ ___ sin. («— b) r' 
P — sin. (tf— *,) ■" T 

cos. y = cos. x cos. x' cos. (c'«— c) -f- sin. x sin. x'. 
sdly. We must reduce into numbers the following formulas, 
leaving p indeterminate, viz. 

f- = R a ■+■ sR cos. x cos. («-t).p-f p a 
r" = R' a -f sR'jQ cos. x' cos. (e'-c*) . P + /3" . / 

V — RR'cos. (e'— <?) + | RjQ cos. x' cos. (<?~c') +R' cos. x 

cos. (^ — c) | . p -}* /? cos « y '• P*» 

3dly. We must determine p by means of these formulas, viz. 
fe 2 = r* + r" — sV 
a 2 = er 2 -f 2/-" 

In the trials necessary for approximating to the value of p, 
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the terms of the final equation, which are of the fourth order, 
being always inconsiderable ; they may be omitted at first, 
and then only taken into the account when a near value of p 
has already been obtained. 

It has already been shown that S = £. ? os ' - ~ » whence £ 

= — -7"s — ;, . cos. x°; therefore, since cos. x° is in every case 

COS. I C "■■■/( 1 v 

positive, the quantity cos r >_ H ) w iH nave the same sign as £: 

and thus, from what is proved in No. 5, we infer that the 
comet's distance from the sun will be less or greater than the 

earth's distance, according as the sign of , , is positive 
or negative; and the two distances from the sun will be equal 

when — T3 — -^ = 0. These observations will often enable us 

cos. (e*— n) 

to assign such first values of p as will lead to a solution with- 
out many trials. 

The preceding method will always give one solution. To 
prove this: let B denote the chord of the earth's orbit drawn 
between the places of the planet at the two extreme observa- 
tions ; and further, let 

R+R'+B 
cos. u = 1 ~- 

2 

, R+R'— B 
cos. u = 1 - : 

2 

then if we apply to the earth's orbit the formula for finding 
the time of describing an arc of an ellipse by means of the 
chord of that arc, and the two radii vectores drawn to its ex- 
tremities ;* we shall get, 

6 = (u — sin. u) — («' — sin. «'). 

Let cos. m = 1 — : then cos. u = cos. m — ~ B ; and 

* Mec. Celeste, Prem. Part. Liv. II. § 27 formule a. 
MDCCCXIV. X 
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cos. u' = cos. m -\- ~ B : we may therefore suppose in series, 

u - sin. u = M< o) - M (I) . 1 + M (z) . £ - M (3) .A" + &c. 

«'- sin. «'=M< 0) + M(°. 1 + M (2) . | + M (3) . | + &c. 

where M*- = m — sin. m; and M , JVT &c. are coefficients 
derived from M*° : hence, 

= M (1) B + M 3) .- + &c; 

4 

and by squaring and neglecting the fourth and higher powers 
of B, we get, 

B* = 



M<0} 2 

but, from the theory of series, we get 



A/1 0) <* • ( m — ' Sln ' *") _ I "~ cos " m , _/ ^' ~ cos - 

~"~ "™" d . cos. m "*~ sin. m "~~ \r l _j_ cos . 



therefore =4- = \ + 2 1' 2 = ( b Y substituting the value of 

cos. m) ^r - i. Let A 2 = aR' + 2 R" = (R + R')' + 

m _ R' )* ; then we may substitute A for R + R'» since quan- 
tities of the fourth and higher orders are neglected in the value 
of B* : consequently we shall have 

But, when quantities of the fourth order are neglected, the 
final equation of the foregoing method will become b % = ^- . 
and B, R, R', A are what b, r,r', a become, when P = o ; thus 

, fix 

it appears that b" in less than —-, when p = o. But as p in- 
creases, at least after a certain limit, 6* will increase with- 
out limit, and ^~ will decrease without limit ; so that, when 
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f = en, b* will be infinitely great, and *£ will be infinitely little. 
Therefore there will always be one value of />, between the 

4.6* 

limits and c/j, that will satisfy the equation 6* = *-. 

10. Before applying the method here proposed to examples, 
it will be necessary to investigate the formulas by which the 
elements of the orbit are to be computed. 

Let u and u' denote the true anomalies of the comet at the 
two extreme observations ; or the angles which r and r' make 
with the perihelion distance of the parabolic orbit ; then in the 
triangle formed by r, r', and b, the angle opposite to b will be 
equal to the sum or difference of u and v': therefore If = r* + 
r* 2 .— arr' cos. (u + u') : but it has already been shown that b* 
==z y- _j_ r n •— sV : therefore V = rr* cos. (u ± u'). Let the 
angle contained by r and b in the same triangle, be denoted 
by <p ; and that contained by r 1 and b, by <p': then it is plain that 

r—r 1 cos. (a±u') , t r'—rcos.(v±v') . 

cos. q> = j-* ; and cos. <p = ^ ■ : whence 

we get the following formulas for finding the angles <p and p', 
which are the angles that r and f make with the chord b, viz. 

r*-V 

cos. <p — - 7 - r 

cos. <p' = - ; . ' ■ : 
r r . b 

observing that regard must be had to the signs of the cosines, 
as well as to their values in numbers. Again let 

— COS. \\) = -7— 

,, r'-r 

— cos. 4/= -y; 

then the angles 4> and 4>', the cosines of which differ only in 
their signs, will be supplements of one another ; and, as it is 
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easy to prove from the nature of the parabola, they will be 
equal to the angles which the chord b makes with the axis of 
the orbit. From the angles thus found, the true anomalies 
are immediately deduced : for u is equal to the difference of <p 
and fy; and u', to the difference of <p' and i}/. Further, the time 
of passing the perihelion will be between the two extreme 
observations, when $ and 4>' are both greater than <p and q>' 
respectively : otherwise the time of passing the perihelion will 
be before, or after, both the observations, according as r, the 
radius vector at the first observation is less, or greater, than r*, 
the radius vector at the third observation. But these rules 
suppose that the angular motion of the comet in its orbit, in 
the interval between the extreme observations, is less than 
180 ; Which, in fact, will comprehend all the cases that can 
occur in applying the method. 

The true anomalies of r and r 1 , being known, the perihelion 
distance, denoted by D, will be found by either of these for- 
mulas, viz. 

D = r cos. 2 — 

D = r* cos. 2 — . 

To find the time of passing the perihelion, we must take 
the times corresponding to v and t/ from a table of the motion 
in a parabola ; these times, being multiplied by D*, will give 
the intervals between the two extreme observations and the 
time of passage. 

In order to determine the position of the orbit in the hea- 
vens, it is best to begin with seeking the heliocentric latitudes : 
let these latitudes be / and /' ; then 

sin. / =b ~ . sin. x ; sin. I' == -j . sin. x'. 
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Let E and E' denote the two angles of elongation (or the 
differences of the geocentric longitudes of the sun and the 
comet) ; and C and C, the two angles of commutation (or the 
differences in longitude of the earth and the comet) : then 

y-, O COS. ?t T, 

sin. C = f . — t x sin. E 

sin.C=»$.g£* xS in.E'; 

by means of which formulas, the heliocentric longitudes of the 
comet at the two extreme observations will be known. 

Let TNn represent the ecliptic; K and K', the places of the 
comet in the ecliptic, determined by the heliocentric longitudes 
already computed; PK and PK', the circles of latitude, and 
KC and K'C the heliocentric latitudes already found ; and lastly, 
let NCC'w be the great circle in which the plane of the orbit 
meets the heavens : then although it is possible that the comet 
may pass from C to C, either by describing the small arc 
CC, or the arc which CC wants of a whole circle ; yet from 
the nature of the case, no ambiguity can hence arise, because 
in finding a first approximation to the orbit of a comet, the 
arc CC, which embraces the whole motion in the orbit in the 
interval between the extreme observations, will never con- 
tain a great number of degrees : it is therefore easy to infer 
whether the comet has increased, or diminished its longitude, 
in passing from C to C ; in the first case, the motion of the 
comet is direct ; in the second case, it is retrograde. Nor can 
there be any difficulty, or ambiguity, in distinguishing the 
Ascending Node, through which the comet passes from the 
south, to the north, side of the ecliptic. The ascending node 
is marked, in the figure, with N ; the descending node, 
with n. 
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The arc CC'is the measure of the angle at the sun's centre, 
contained by r and r'; and it is known by means of the true 
anomalies, since it is = v + v' ; or, by means of the formula, 

cos. CC = jp: PC and PC are the complements of the he- 
liocentric latitudes; and KK/ is the difference of the heliocentric 
longitudes. Now in the spherical triangle CPC, we have 

sin. CC : sin. P, or sin. KK' : : sin. PC, or cos. /' : sin. C ; 
and in the right-angled triangle KCN, of which the angle N is 
the inclination to the ecliptic, denoted by I ; we have 

Rad. : cos. KC, or cos. / : : sin. C : cos. I ; 
therefore by combining these two proportions, we get 

T sin. KK' 7 7/ 

COS. I = -. — r^r, . COS. / COS. / . 
sin. CL. 

The resolution of the right-angled spherical triangle KCN, 
will give the longitude of the ascending node, and the place 
of the perihelion on the orbit. For the side KN is the differ- 
ence in longitude between the comet at the first observation 
and the ascending node: and the hypothenuse NC is the an- 
gular distance on the orbit between the same node and the 
comet at the first observation ; and from this it is easy to find 
the angular distance between the node and the perihelion, 
which will fix the place of the perihelion on the orbit. 

11. In applying the preceding method, the examples in 
Legendre's Memoir and the supplement to it, have been pur- 
posely taken : first, because these make a selection of instances 
greatly varied in their circumstances : secondly, because the 
results obtained by the formulas investigated in this paper can 
thus be compared with his, and other methods, and with the 
corrected elements, which are likewise given in his Memoir. 



a first Approximation to the Orbit of a Comet. igg 



Application to the second Comet of 1781. 

The following observations of this comet, which was dis- 
covered by M. Mechain, are taken by Legendre from the 
Mem. de I' Acad, des Sciences, for 1780:* they are all reduced 
to the same hour of every day, namely to 8 h 29' 44" mean time 
at Paris. 



Times of 
Observation. 


Longitude. 


North Latitude. 


Longitude of ^. 


Log. of R. 


September. 


c, 3°7 H 45 


*, 55 17 9 


e, 52 54 « 


R > 9'994 86 4 


19 . 


c°$o6 51 26 


*»° 39 14 48 


e>°57 57 4R. 9994426 


24 


c/306 42 20 


x,' 31 4 52 


£,'63 41 


R/ 9'994°28 



By applying the directions and formulas of No. 8, 1 have 
found, 

Log. 46*= 0/0732228' 
n = 306°- 19' 5" 

i =as 89 21 27 

£ = — o- 0000970 

» = 39° 1 4' 49" 
A = 5£ 17 28 

/*'= 31 5 o 

log. /3 = 0*2892147 

cos. <y = 99600258 

With regard to the three principal formulas of the orbit, I 

remark that it is the logarithms of the coefficients, and not 

the coefficients themselves, that are afterwards used : on this 

* Nouvelles Methodes, p. 33» h x * x ii- 
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account it is convenient to denote the coefficients by means of 
their logarithms : which is done, by writing, for instance, 
num. (9.4,824,80) for the number whose logarithm is 9-482480. 
The formulas of this orbit, expressed in the manner here pro- 
posed, are as follows, viz. 
r* = 0-976625 ■— num. (9*482480) xp + p' 
r' 2 = 0-972873 — num. (0-163532) x p + num. (0-578429) x p B 
V = 0-959610 — num. (9-847341) x p + num « (0-249241 )xp 2 . 
We are now prepared to approximate to the value of p; for 
which purpose we must seek a first value as near the truth as 

we can. Now, in this instance, ^7^5^ is positive, and very 

small : hence, we may infer that r will be less than R, but, 
approaching near to an equality with it. Let us then suppose 
f- = R 2 = 0*976625; and we get p = num. (9-482480) as 
0*3037 ; or p = 0-3 nearly. By substituting 0-3 for p, we get 

r 2 = 0-975508 

^=0-876635 

V = 0-908288 

b*=xr* + r'* — 2V — 0-035567 

a* = ar* 4- ar" = 3-704286 . . log. 0*5687044 

log. a = 0-2843522 

■If — € — __ 0-025932. 

The error here having the same sign as in the supposition 
p =a (No. 9) ; we hence learn that p is greater than 0*3. 
Suppose p = 0-4, then 

r* = 1*015136 

r' l = 0*996083 

V= 0962188 
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b' = 0*086843 

a* = 4-022438 . . . log. 06044892 

log. a = 0-3022446 
f— £ = + 0-027826. 

This error is nearly equal to the former, but it has an op- 
posite sign : hence p will be nearly = 0-35. Suppose p = 0*35, 
then 

r 4 = 0992821 

r'*= 092688 8 

V = 0-930800 

£>•=: OO58109 

a" as 3-839418 . . . log. 0-5842654 

log. a = 02921327 

6 2 — *j- = — 0-002298. 

By comparing this error with the last, we get a nearer value 
of p, viz. p = 0*354 : and the approximation is now far enough 
advanced to take in all the terms of the final equation in making 
the next substitution, viz. p = 0-354 : then 

r» = 0994422 

r- l = 0-931729 

V = 0-932985 
b* = 0-060181 

a' = 3-852302 . . . log. 0-5857202 

log. a = 0-2928601 
1 + i • ( — ] + Ti • ? = 1-001434 

5» _» 2_ x t -001434 = — O-OO0212. 
MDCCCXIV. Y 
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In order to get another equation to compare with this last, I 
resume the equation of the supposition, p = 0*550, and correct 
it by taking in all the terms of the final equation : then 

b* — — x 1-001434 = — 0002385 : 

and now by comparing the two errors I finally get p = 0-35439. 
This value being substituted, the result is as follows, viz. 

r* = 0-994580 

r'*= 0932205 

V = 0-933200 

6 2 = 0-060385 

<*" = 3 -8 5357° • • • log- 0-5858633 

log. a = 02929316 

4.6* 
b* — *- x 1-001434 == + 0-000002. 

Having now obtained the exact values of p, r", r /l , V and ff, 
it remains to deduce the elements of the orbit from them. 
Now we have (No. 10) 

£~ == cos. <p v <p = 75 29' 43" 



r—-f 



= — cos. xj/ -.- xj/ ===== 97 25 50 



Anomaly u == 21 56 7 

r**—V 

-py = cos. <p' v <p'= 90* 14' 25" 
= — cos. 4>' '.• 4/= 82 34 10 



6 



Anomaly u'== 7 40 15. 
For the perihelion distance, we have 

D = r cos.* ~- •.• log. D = 99828083 

D = r'cos.* — -.• log. D =5 9-9828085 
D = 0-961188. 



a first Approximation to the Orbit of a Comet. 16% 

By a table of the motion in a parabola I have found the in- 
tervals corresponding to u and v', as follows, viz. 

16-1316 days, and 5-5198 days 
log. 1-2076773 07419233 

log. D* = 90742 123 9-9742123 



0-7161356 
Num. 5-2016. 



1-1818896 
Num. 152016 

These are the intervals between the two extreme observa- 
tions and the passage of the perihelion, which according to 
the rule in No. 10, will be posterior to both the observations: 
therefore the comet will be in the perihelion, Nov. 29' 5556, or 
Nov. 29, 13 14 20' 4." 

M. Mechain has applied the method of Laplace to this 
comet, using as the basis of his calculations, five observations 
between the 24th and 25th of November.* Laplace's method 
gives rules for finding two elements only, viz, the perihelion 
distance, and the time of arriving at the perihelion : for the 
sake of comparison, I now subjoin the results of M. Mechain ; 
the same elements calculated by Legendre from the obser- 
vations used here ; and the more exact elements corrected by 
distant observations. 





Perihelion 
Distance, 


Time of Pas- 
sage Nov. 29. 


M. Mechain by Laplace's method 


°'95H5C>9 


h 

18 10 34 


Legendre 


0-9604,49 


16 18 19 


In this paper - 


0961 188 


13 20 4 


Corrected elements 


°-9 6o 995 


12 42 46 



These calculations are given in Biot's Astron. 2d Edit. Vol. III. 
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I have likewise computed the remaining elements of this 
orbit by the rules laid down in No. 10, and I have found them 
to be as follows, viz. 



Inclination 



2 7° i' oo" 



Place of the ascending node J? 56 32 
Place of the perihelion - 15 59 18 
The method here proposed will represent the two extreme 
observations exactly ; so that all the error of the solution will 
fall upon the middle observation. I have therefore employed 
the elements found above, to compute the geocentric place of 
the comet, Nov. 19, 8 h 29' 44,", in order to compare it with 
observation. The result of this calculation is as follows, viz. 



Calculated 
Longitude. 



306 51 30 



Observed 
Longitude. 



306 51 26 



Error. 



+ 4 



Calculated 
Latitude. 



Observed 
Latitude. 



39 *4 3439 14 4 8 



Error, 



•H 



In Legendre's memoir I find a like comparison of the cor- 
rected elements with observation; the error in longitude being 
+ 100"; and the error in latitude, — 55". It appears then 
that the errors of the elements found here are almost nothing, 
and even much less than those of the corrected elements. 

Application to the Comet of 1 769. 

The places of this comet set down below are not derived 
from actual observation; they are calculated by Legendre 
from the known elements of the orbit.* In this instance there- 
fore, the inaccuracy of the results is to be attributed, not to 



* Nouvelles Methodes, p. 43, § XL. 
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the errors of observation, but solely to the defects of the 
method. 



Times 
Sept. 


Longitudes. 


South Latitudes 


Longitudes of®. 


Log. R. 


d. h. 

8 14 


• 1 n 

c, 101 18 8 


x, 22 14 35 


e 34 6 35 31 


R, 0-0026648 


IO 14 


c°, 1-12 51 31 


x°,23 28 15 


^,348 32 22 




12 14 


c', 124 26 47 


a',23 4 8 3 6 


^350 29 20 


R', 0002 1838 



From tjhese data, I have found 

log. 4# 2 = 8-2773428 
n = 33 22' 56" 
i == 23 48 49 

I = -{- 0-0006684 

« = 8o° 23' 45" 
k = 6g 38 2i£ 
#=90 58 25 
log. |Q':= 0-0071924 

COS. <y = 99691708. 
1* = 1*012347 — num. (9-8913194) x p -f p* 
r"= 1-010107 — num. (0-1132123) xp + num. (0.0143848) xp* 
V= 1-008888 — num. (0-0104917) x p + num. (9*9763632) x/>*. 
In this instance — ~ — ? is positive : therefore the comet is 

COS. («• — H) l 

nearer the sun than the earth : hence p z. num. (9' 8 9 1 3 1 94)» 
or p 0778. In order to get narrower limits, I take p = 
= o*3 8 9> or p = 0-4, for a first approximation ; then 



0-778 
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t* == 0860903 

t* = 0-656370 

V = 0.750632 

f = r* + r' % — 2V = 0-016009 

a" = 2r* + ar" = 3*03454,6 . . . log. 0-4820936 

log. a = 0-2410468 
& — — = + 0-60513?. 

As this error is positive p must be diminished: therefore 
suppose p = 0-3 ; then 

r' = 0-868764 
r"=z 0713794 

V = 0-786784 
6* = 0-008990 

a" = 3-165116 . . . log. 0-5003896 

log. a = 0-2501948 

fel — IT = "~ °* ool6 55- 
By comparing the two errors it will appear that p is between 
0-32 and 0-33: suppose p = 0-33, and because the errors are 
now small, let all the terms of the final equation be taken into 
account; then 

r 2 = 0-864306 
r" = 0-694396 

V = 0-773950 
b* = 0*010802 

a* = 3-117404 . . . log. 0-4937930 

log. a = o- 2468965 
1 +'*(—) + » * ? = rool 774 
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&■■—£- x 1-001774, = -f- 0-000057. 

Because the last value of p is too great, let p = 0-329 : then 
r* = 08644,26 
/•== 0-695012 

v = ©77435 1 
b* = 0-010736 

a* = 3-118876 . . . log. 0-4939875 

log. a = 0-2469938 
6« _-. 1_ x 1*001774 = — 0-000007. 

By comparing the two last errors we finally get p =r 
0-32911 : then 

r* = 0-86441 2 
r"= 0-694945 
V= 0774307 
6*= 0-010743 
<f = 3-118714 . . . log. 0-4939760 

log. a = 0-2469880 
If— ^ x 1-001774 == 0. 

From the values just found, we get 

log. r — 9-9683604 

log. r* =9-9209752: 
and the angle contained by r and r ', the cosine of which is = 

,-jp, will thence be found = 2 31' 35''. These quantities are 
sufficient for determining the parabola described by the comet ; 
and as they are the immediate results of the method, the fairest 
way of judging of the exactness of that method, seems to be, 
to compare them with the true values calculated from the 
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known elements of the orbit. Now the elements employed 
by Legendre In computing the places of the comet used in 
this example, are these, viz. 

log. of perihelion dist. g 030847 
Time of passage, Oct. 7*5310; 

from which data I have calculated the radii vectores and the 
true anomalies for the 8th and 1 2th of September at i4 h , as 
follows, viz. 



September. 


Log, of Rad. Vect. 


Anomaly. 


d. h. 

8 14 


9968361 


137 17 35 


12 14 


9.920961 


134 46 10 


Angle between r and r' 


2 31 25 



These quantities differ very little from the results obtained 
above : the errors are indeed hardly greater than the dis- 
crepancies which it is difficult to avoid in a complicated cal- 
culation, viz. first computing the geocentric longitudes and 
latitudes, and then going back from these to the elements of 
the orbit. 

The exactness of the approximation is sufficiently proved 
by the comparison already made : but I have likewise calcu- 
lated these two elements, viz. 

log. of perihelion dist. - 9-091920 

Time of passage, Oct. - 7*5403 

the perihelion distance being 4^th part too much, and the 

time of passage about 13^' later than the true time. These 

errors, although very .small, are yet greater than might be 
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expected, considering the exactness with which the comet's 
distances from the sun, and the included angle, were deter- 
mined : but a very little variation of these quantities will pro- 
duce a great alteration in the magnitude of the parabola, when 
the angle contained by the radii vectores is so small as it is in 
this instance. 

Of the two comets to which the method in this Paper has 
been applied, the first varies very slowly in longitude, and 
very rapidly in latitude ; the second, on the contrary, has a 
quick motion in longitude, and a very slow motion in latitude : 
in both instances however the intervals between the observa- 
tions are equal ; and I now proceed to give other examples 
where this condition is not observed. 

Application to the second Comet of 1805. 
The following observations are the same that Legendre 
has used for determining the orbit of this comet :* they were 
communicated to him by M. Bouvard. 



Times 
November 


Longitudes. 


North Latitudes. 


Longitude of ©. 


Log. R. 


23-32241 


• 1 a 

c, 24 41 5 


x, 27 25 19 


a > 11 

e, 61 8 23 


R> 9'994P°4P 


30-51095 


c°,i5 3 8 3 6 


\°,ig 25 6 


e°,6S 25 47 


R°>9 993 66 55 


December 
5'2958l 


c', 2 6 33 


a', 3 19 22 


e',73 17 5 


R', 9-9933784 



Calculating from these data, we get 
log. 46' = 92296576 
n= 359° si' 58* 



at 



* Supplement aux nouvelles Methodes, $ 30, p. 36. 
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i= 50 30' 11" 

§ = — 0*0039318 

a — 25 21' 33" 

k = 36 38 32f 

^'= 4 18 27 
log. t = 0-8566407 
log. •/= 0-6798693 

log. /3 = 9-5593 8o 3 
cos. 7 = 9-9268248 
r* = 0*973662 + num. (0-1489016) xp + f 
/•= o- 96gg6j + num. (9-3618159) x p + num. (9*1 187606) xp a 
V = 0-950062 + num. (988 19026) x p + num. (9-4862051 )xp 2 . 
In this instance, the form alone of r 2 shows that r > R ; and 
likewise that r increases continually as p increases : we will 
therefore make different substitutions for p, increasing the 
value when the error is negative, and decreasing it when the 
error is positive. 

First, let p = o-i : then 

r 2 = i- 124556 

r"= 0-994286 

V= 1-029316 

b* = r* -j- r* — 2V = 0-060210 

a* = 2r* -}- sr' 2 = 4-237684 . . . log. 0*6271286 

log. a = 0*3135643 
*•_£ = _ 

a 

Next, let p = 0-2 ; then 

r *= 1-295456 
r''=z 1-021234 
V == l- 1:14698 
b % = 0*087294 



a 0-022221. 
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a* = 4-633380 . . . log, 0-6658978 

log. a = 0-3329489 
lr — € = + 0-008461. 

A comparison of the two errors will show that p = 0-17 nearly: 
this value being substituted, we get 

r 2 = 1-242087 

r' 2 = 1-012874 

V= 1-088439 

b % = 0-078083 

a* = 4-509922 . . . log. 0-6541690 

log. a = 0-3270845 
b 1 — — = — 0-001822. 

a 

By comparing this error with the last, it will be found that |» 
is between 0-175 and 0-176 : I therefore put p = 0*176; and 
as the errors are now small, I employ the complete final equa- 
tion: then 

f = 1-252617 

r" = 1-014527 
V = 1-093647 
If = 0-079850 

a* = 4 -534 s88 • • • lo S- o6 5 6 509i 

log. a = 0*3282546 

1 + I • \(~r)' +£•■£= 1-002846 

b* — — x 1*002846 sss — 0-000067. 

I next put p = 0*177 : tnen 

r* = 1-25437 8 
r' 2 =-= 1-014803 

Zs 
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V = 1-094516 
b*= 00801 49 
•**= 4*53 8 3 6s • • • Io g- 0-6568991 

log. a = 0-3284496 

h " "~ T x 1,0028 4 6 — + 0-000268. 
From the two last errors we finally get p = 0*17620 ; which 
must be substituted : then 

r a = 1-252969 
r' 2 = 1 -01458a 
V= 1-093821 

If S3 0-079909 

« 2 = 4-535102 . . . log. 0-6565888 

log. <z = 0-3282944 
6 s __. i, x 1-002846 = 0. 

From these values we get 
~~ = cos. <?> •.• <p = 59° 48' 14" 

r"- y 

-pf- = cos. <p' -.- <p' = 106 9 30 

r-r* , , ., f-i/ = 119° ai' Aft" 

~ = -<-•* = + -, 4/ v{* /== J 38 ^ 

anomalies { u = * ~ * = ^S' 3*'' 
U = <?>' — ip = 39 31 15 

Then, to find the perihelion distance D, we have 
D = r cos.* ~ = 9-9504285 logarithms 

D = /COS.* ~ as 99504285. 

The times corresponding to the anomalies, taken from a 
table of the comets, are as follows, viz. 
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days 45"Oi35 days 30-8044 

log. 1-65334,28 log. 1 -48861 27 

log. Dl 99256427 log. Dl 9-9256427 



*'57 8 9$55 1-4142554, 

Num. 37-9303 Num. 25-9570 

hence the time of passing the perihelion is fixed for Dec. 
31-2527. 

By making the rest of the computations in No. 10, 1 have 
found the following elements of this comet, viz. 
Log. perih. dist. - 9'95°&9 
Time of passage, Dec. 25-2527 
Inclination - i6°3o' 37^-" 

Place of ascend, node 250 33 34 
Place of perihelion 109 20 21 

Motion Direct. 

Legendre, from a comparison of all the observations of this 
comet, which are only six in number, has found the following 
corrected elements, viz. 

Log. perih. dist. - 99502700 
Time of passage, Dec. 25-28551 
Inclination - i6°3i' io'' 

Place of ascend, node 250 33 34 
Place of perihelion 109 23 39 
which differ very little from those found above, and prove the 
exactness of the approximation by the method in this paper. 

Neither of these two systems of elements represent the ob- 
servation of Nov. 30, with much accuracy, the errors in lon- 
gitude amounting to between -j~ a'-|- and -f- 3' ; on which 
account it is probable that there is some peculiar inexactness 
in that observation. 
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Application to the first Comet of 1805. 

This comet was discovered by M. Bouvard : and Legendre 
has given us five observations communicated to him by that 
astronomer.* Of these observations those which I have se- 
lected for the purpose of computing the orbit, are unfavour- 
able in several respects. For, besides that the intervals between 
them are too unequal, the whole time elapsed, which embraces 
a heliocentric motion of no less than 25 , is rather too long : 
and it happens that the comet is placed with regard to the 
earth and the sun, so as to approach very near the limit, when, 
as has been remarked in No. 9, the method will fail. A better 
choice cannot however be made from the observations recorded 
by Legendre, without interpolating. 



Times. 


Longitudes. 


North Latitudes. 


Longitudes of©. 


Log. R. 


October. 
22-6849 


C, 163 20 53 


*» 22 59 53 


e, 29 19 5° 


R » 9'9974>* 1 


30-6867 


^,183 48 32 


x°,i5 37 21 


e °>37 39 4# 


R°> 9*996496 


November. 
37209 


c', 191 46 15 


\',i2 2 29 


e',41 22 26 


R', 9*996053 



In this instance, we get 

log. 40* == 9*2341873 

rc = 34 58'3ii" 
1 = 28 25 49 
$ = + 0-0246945 
u =s 146 17' 25" 
h == 124 51 1 
h' = 154 o 41 

* Supplement, p. 14, § 15. 
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log. t = 0-9031887 
log. / = 0-6057574 
log. /3 = 0-1372458 

COS. y = 9-94,11400 

y* = 0-988192 — num. (0*104392) . p + p* 

r* % = 0*981987 — num. (0-363921 ) . p -f- num. (0-274492 ) . p* 

V = 0-963404 — num. (0-244138) . p + num. (0-078386) . p 2 . 

In this instance — ;-= — r is positive ; and therefore the comet 

cos. (e°— n) r ' 

is nearer the sun than the earth is : hence p < num. (o- 104392), 

or p < 1.270. To get narrower limits I assume p = -~~ = 

0635, or p = o-6 ; then 

r 2 = 0*585160 
r^ss 0*272317 
V == 0*341949 
fc 2 — 0-173579 
a 2 = 1714954 • • • lo g* 0-2342388 

log. a = 0*1171194 

4.9* 
6* — ±- = -j_ 0*042642. 

This error being positive, the value of p must be diminished : 
therefore suppose p = 0*5 ; then 
r* = 0-602332 



r"= 


296528 






V== 


03^5 G 35 






6 2 = 


0-127590 






«'== 


1*797720 . . . 


. log. o-< 


2547220 




log. 


a= o-i 


ia/3610 


*>*- 


4.6* 

- 2- =-= — OOO297. 
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Here the error being very small, it will be proper to take in 
all the terms of the final equation : then, 

1 + £•(—) x n . -5=1-080389 

4.6* 

b* — — x 1-020382 = — 0-002903. 

It now appears that p is between o-6 and 0-5, but much nearer 
to 0-5 : therefore put p = 0-51; then 

r*~ °'5997H 
r /2 = 0-292413 

V = 0*380190 

h*= 0-131747 

a 2 = 1784254 . . . log. 0-2514566 



log. a = 0*1257283 

/r x — r'*\ 1 b z 

1 + i * (-T-J + Ti • T> = *'°*°9H> 

b 2 — ^ x i'20984 = -f- 0.000685. 

By comparing this error with the last we get p = 0-5081, which 
must be substituted ; then 

r 2 = 0*600197 

r»= 0*293163 

V = 0-381205 

b 2 =z 0*130950 

a" = 1*786720 . . . log. 0*2520565 

log. a = 0*1260283 

4.9* 
fr 2 = -^ x 1*02086 = — 0-000005. 

The approximation has already been pushed farther than 
was requisite : for, in this instance, the terms of the fourth 
order affect the errors in the third place of decimal figures ; 
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and therefore the terms of the higher orders, which are omitted, 
may be supposed to affect the errors in the fourth and fifth 
places of figures, beyond which degree of exactness the ap- 
proximation is already carried. To Compute the elements of 
the orbit, we have 

£=%■ = cos. *vf = 38*38' 6" 
r ~ = — cos. 4> v ip = 130 8 21-§ 



Anomaly v — $ — <j> = 91 30 15! 
^=f = cos. <?>' -.- <p' = 1 16 42' 6£" 

~ — _ cos. 4/ v 4;' = 49 51 38*- 
Anomaly 1/ = <?/ — 4/ = 65 ,50 28. 
Then, D = r cos. 2 — = 9-5765635 logarithms 

D == r' cos. a — = 9-5765635. 

The times corresponding to the anomalies, taken from the 
table, are 

days 114-04,8 days 62*1248 

log. 2-0570877 log. 1-7932650 

log. D 1 , 9-3648453 log. D 1 , 9'3 f5 4 J 8453 

1 *4* 1933° 1-1581103 

Num. 26-4200 Num. 14-3916: 

the difference of these two intervals is 12 0284 ; but it should 
be 12-0360: I therefore add 0-0038 to the first, and subtract it 
from the second interval ; then they are 26-4238 and 14-3878: 
and by adding these corrected intervals, to the times of the 
observations, the time of the comet's passing the perihelion 
will come out, Nov. 18-1087. 
mdcccxiv. A a 
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Having completed the rest of the calculations in No. 10, the 
elements of this comet will be found as follows : viz. 
Log. of perih. dist. -?. 9 576563$ 
Time of passage, Nov. - i8'io87 
Inclination - - 15° 40' 17" 

Place of ascend, node - 34,5 7 30, 
Place of peri. - 148 11 7 

Motion direct. 
These elements will not represent the two extreme observa- 
tions with perfect accuracy as in the other examples, on account 
of the correction which it was necessary to introduce, in order 
to annihilate the difference between the time in the orbit and 
the actual interval between the observations : but the errors 
proceeding from this cause will be inconsiderable, and can 
hardly amount to 20." If we employ the same elements to 
calculate the place of the comet, Oct. 306867, we shall find 

these results, Viz. Geo. Long. Geo. Lat. 



Calculated i83°5i'20" 
Observed 183 48 32 



i5°3 6 '37" 
15 37 21 

error — 44 



error + 2 48 

These errors are in reality less than might be expected, 
considering the disadvantages under which the computation of 
this orbit is made. Legendre has given a system of elements 
which represents the observations of the comet very exactly, 
the greatest error not exceeding 1': and, for the sake of com- 
parison, I likewise subjoin these elements, viz. 

Log. of peri. dist. - 9'5747<)8 

Time of passage, Nov. - 18-01736 
Inclination - - 15° /,8' 12" 
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Place of ascend, node - 345 6' 51" 
Place of the peri. - 149 o 28 

12. The method of finding a first approximation to a comet's 
orbit, investigated in this Paper, has now been applied to four 
examples in very different circumstances from one another ; 
and the results obtained in them all are satisfactory. These 
instances will suffice for making the rules of calculation clear, 
and for confirming the accuracy of the analysis from which 
they have been deduced. It remains that we now consider 
the particular cases, already noticed (No. 7), in which this 
method fails. 

If a new planet should be discovered that moved exactly in 
the plane of the ecliptic, three geocentric observations of it 
would not be sufficient for finding its orbit. The latitudes 
being wanting in the circumstances here supposed, the longi- 
tudes alone would not furnish conditions enow for determin- 
ing the magnitude and position of the curve which the planet 
described. The same thing will likewise happen in another 
situation, when the latitudes, although not evanescent, yet 
depend, all of them in the same way, upon the longitudes ; 
and that is when the three geocentric places of the planet are 
situated in one great circle of the heavens which cuts the 
ecliptic in the points occupied by the earth and the sun at the 
time of the middle observation. This last case may indeed be 
considered as including the former one, when the plane of 
the planet's orbit coincides with the ecliptic. 

With regard to the comets it must be remembered that three 
complete observations are more than sufficient for finding an 
orbit j and on this account there are still conditions enow, but 
no more than enow, to determine the elements sought, in the 

A a 3 
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circumstances abovementioned in which the problem becomes 
indeterminate when applied to the planets. We shall therefore 
proceed to analyze this case, which has already been noticed as 
forming a real distinction in the problem we are considering. 
Resume the three following equations which have already 
been investigated in No. 4,, viz. 
<r'p cos. A -f- <r/ cos. ti — <r"f cos. x° cos. (c — n) = 

Tt'(t+t') / I I \ TJQ f O \ 

■ T- • [FS - rs) • R cos. («°— ») 

e-'p cos. i sin. A -{- c?' cos. z sin. A' — <r"p" cos. x° sin. (c°— n) = 

-=^ 2 .(i-E 3 ).R°sin.(^-i.) 

<r'p sin. / sin. A, -f" *?' sm - z sm « ^' — ^V sm< *°= : 
and employing A° and i° to denote the same things as in 
No. 5, let the following values found there, viz. 
cos. x° sin. (c°— n) = cos. i° sin. A° 
sin. x° = sin. i° sin. A° 
be substituted in the two last of the equations above ; then 
these equations will become, 
v'p cos. i sin. A -J- <rp' cos. i sin. A — <r"p° cos. i° sin. h° =■ 

4— ' (p — 53] • R sin - ( e — «) 

er'p sin. z sin. A -j- 07/ sin. z* sin. A — <r"p° sin. i° sin. A° = : , 
further let these two equations be added together, after being 
multiplied by cos. i and sin. i respectively : then we shall get 
c-'p sin. A + <rp' sin. A' — <r"p° cos. (z — z'°) sin. A = 

- '"—^ • (£ ~ &) R° cos. i sin. (*•-*), 

In the particular cases we are considering, either i = z°; 

or i and z° are both evanescent ; in both which cases cos. (f— 4°) 

= 1. And if we attend to the analysis in No. 5, it will appear 
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that, in every case, sin. {i — z'°) is of the second order; and 
hence cos. (z' — z* 9 ), differing from 1 only by quantities of the 
fourth order; may generally be regarded as equal to unit. The 
last equation will therefore become, 
<r'p sin. h -{- <rp' sin. H — <r''p° sin. h° = 

tt'(t+t') /I I 1 n . • • / o v 

~T~ • [75 — w) • R cos - l sin ' ( fe — n ) : 

and, if this equation be combined with the first of the three 
equations set down at the beginning of this article,* so as to 
exterminate p° and p', we shall get, 

cr'psin. {W-h) - <r P ' sin. {K-h°) = - =^2 . (J, - J- 3 ) .R° x 
{ cos. (e° — n ) sin. A° — cos. i sin. (<?° — n) cos. /z°| ; 

*> sin. (/*'-/*) _^°sin.(A'-/i ) = -I^±l2. (£ -^).R°x 
J cos. (e° — n) sin. A' — eos.z sin. (<?° — «) cos. h'\» 

Observing then that — = — , and -^ = — 7 , nearly; if we put 

r> y(r+y) rio cos. (e°— a) sin. £i°— cos, t si n. (V— w) cos, b" ^ 

2 ' ' sin. b'—b 9 ' ~~~ '' 

_ sin, (b ' —b) t' t 
* — sin. (b'—b*) ' 7+P » 

we shall get, 

p = e . p 

,_ sin. (&**-;5) _£_ , JB_ _B 

P — P X sin. (yb'-ib») * t ' r>» R°s : 

the quantities of the second order being omitted in the expres- 
sion of p°; because p° is only to be used in valuing — which is 
already of the second order. 

If we substitute the values of the coordinates (No. 3), in the 
expressions r° 2 = x" 2 -j- y"* + z° 2 ; r 2 = x 2 -f- J 2 + s a ; v' % =- 
x' 2 + y 8 + z' 2 ; V = .rx' 4- jy/ -j- zz'; we shall get 

* In this equation cos. b* must be written for its equal cos. *° cos, (e'— «),- 
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r** = R°* -f. eR° cos. X° cos. (e°— c°) . $ . p + 6 a , P a 
r* = R 2 -f 2R cos. x cos. {e — c) . p + p 2 
/ 2 = R' 2 + 2 R' cos. V cos. («'— c') . / + ? ' a 
V = RR' cos. (<?'— + R cos. x' cos. (<? — <?'). g' 

-{- R' COS. X COS. (<?' — £).£-{- COS. y^'. 

Cos. y denoting here the same thing as before. The first of 
these formulas, which determines r° 2 when o is given, will 
enable us to compute g>' when the value of is known or as- 
sumed, by means of this formula found above, viz. 

, sin. (A 9 — A) f' , B B 

? — ? X sin. (b'—b°) * r •" r° s Ro3' 

Thus the values of the functions, r 2 , r' 2 and V will depend only 
upon one unknown quantity, namely g> ; which may therefore 
be found by the help of the same final equation as in the for- 
mer method. 

The preceding analysis leads us to the following method for 
determining the orbit of a comet, viz. 

1st. We must begin with computing the values of &,n, and 
/, as in the former method. We must then calculate 
cos. h = cos. x cos. [c — it ) 
cos. h° = cos. x° cos. (c° — n) 
cos. h' = cos. x' cos. (c' — n) 

. ; sin. (A'— A) t 

S sin. (A'— b°) ' Z+P 

R •— ll!lil2 zos ' ( £ °— ^) si"- &° — cos - ' sin. (« 6 — n) cos. 2> a p S 
2 * sin. (A'— A 8 ) * K 

COS. y = COS. X COS. X' COS. (</ — c) -f- sin « x SHI. X'. 

sdly. We must reduce into numbers the several coefficients 
of the formulas, viz. 
r» 2 = R° 2 -f 2R cos. X° cos. (e°— c 9 ) . 1 x ? + £ l x f 

, sin. (A — A) _£ , B_ J3_ 

f ~™ ? X sln.(b'-b') ' t • r o3 """ Ro3 
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r* = R 2 + Q R cos - x c <> s ' ( e — ' • f + § 2 

^= R' 2 + sR' cos. x' cos. {e' — c') , ? ' + ? /a 

V = RR' cos, (e'— <?) -j- R' cos. x cos, {e f — c) . ^ 

+ R cos. X' cos. (e — c') , (f 4- cos. <y . ^'. 
gdly. We must employ the same final equation as in the 
last method, to determine ^ , upon which the values of r% r*% 
and V all depend. 

Application to the first Comet of 1805. 

If we turn to the calculation of the orbit of this comet, by 
the former method it will appear that the intersection of the 
great circle which passes through the two extreme places of 
the comet with the ecliptic is very near the place of the earth 
at the middle observation ; for the distance of these two points* 
or the arc c°— n, amounts to no more than 2 21' 16"%. This 
instance therefore is very near the limit when the former 
solution becomes indeterminate ; which takes place when the 
arc e° — n is evanescent, or when the points abovementioned 
coincide. In these circumstances we can hardly ex pect that 
the orbit will be determined with much certainty : for, on ac- 
count of the small divisor tan. [e°— - n), the values of <?, and 
the angle &>, will suffer considerable variations when small 
changes are made in the observed places of the comet, more 
particularly in the mean place. This comet therefore furnishes 
a very proper example for applying the method we have just 
investigated. 

Using the same observations as in the former calculation, 
we shall get, 

log. #=0/2341873 

»= 34°58'3i"i 
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z== 28 25 49 

h = 124 51 1 
h°— 145 29 34 
A 7 =154 ° 4 1 

lOg. B = 0-042397 

log.fi = 8-457324. 
cos. 7 = 9-941140 
And the principal formulas of the orbit will be 
r° 2 = 0-98399 — num. (0-244576) . p + num. (0-084795) . ( 
p' = num. (0-079127) ? + num. (8-457324) . ± — 0-029365 
r 2 = 0-988192 — num. (0-104392) . p + (? 
p n = 0981987 — num. (0-226675) . p' -f- p' 2 
V = 0-963404 — num. (9-683982) . g 

— num. (9-967035) . p '+ num. (9-941140) f ^ 
In order to shorten calculation, I shall assume p = 0-51 
which is very near the true value, as will appear from the 
former computation of this orbit: then 

r°*= 0*40450 . . . log. 9-6069185 

9-8034592 



log. r° 3 = 9-4103777 
log. e '= 9841343 
r*= 0-599714 
r H = 0-294039 
V = 0-382873 
b* a- r * -f- r n -~ 2V = 0-128005 
a" = gr 2 -f- %)■'* = 1-787506 . . . log. 02522476 



log. a = 0-1261238 

, + JL . (~^~j*-f- ± . b ~ = V02O59 

£ a — — x i'Q2059 = — 0-002888. 
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I next put ? = 0-52 ; then 

r° 2 = 0-39946 ...log. 


9-6014733 
9-80O7366 


log. r 03 == 
log. (f = 9- 85008 8 
r 2 = 0-597298 


1 9-4022099 


r' 2 == 0-290041 




V = 0-377429 

b 1 = 0-132481 

a*= 1774678 .. 


. log. 0-2491196 


log 


;. a = 0-1245598 

-I- = 1-02121 


J 2 — t- x 1-02121 = 

(2 


-f- O.OOIO37. 



By comparing the two errors, we shall get £ = 0-5174; which 
must be substituted ; then 

T* = 0-40075 . . . log. 9 6028735 

9-8014367 

log. r^ = 9.4043102 

Io g- ?' = y847 8 3i 
r 2 == 0-597907 
r' 2 = 0-291041 
V = 0-378823 
6 2 = 0-131304 
a 2 = 1-777896 . . . log. 0*2498820 

log. a— 0*1249410. 
If we interpolate 1-02059 and 1-02121, which correspond 
mdcccxiv. B b 
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to the two preceding values of g>, we shall get 1*02106 for the 
present value of the same quantity; then 

b* — — x 1 '02106* =s 0. 

If we calculate the orbit from the values of £, %, f 2 , r' 2 and b" 
that have just been found, we shall get the following ele- 
ments, viz. 

Log. of perihelion dist. - 9*574910 
Time of passage, - Nov. 18 0262 
Inclination - - 15°55'47" 

Place of the ascending node 345 2 8 
Place of the perihelion - 14855 7 
These quantities-come surprisingly near the corrected ele- 
ments of Legendre : and, without further calculation, we may 
presume that they will represent even the middle place of the 
comet, Oct. 30th, within a minute of a degree. 

Although the method of this article has been investigated 
for the purpose of supplying the deficiency of the former 
method, it is nevertheless quite general, and will apply in all 
cases whatever: and perhaps of all the methods hitherto pro- 
posed, this is the one in which the results will be least affected 
with the errors of observation. Perhaps too we are warranted 
in drawing this conclusion from the contents of the preceding 
paper, viz. that a first approximation to the orbit of a comet 
may be deduced from three geocentric observations that shall 
come nearer the true elements than has usually been thought. 
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